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SUMMARY 


The theory for single-degree torsional instability of a two- 
dimensional wing oscillating in a supersonic stream is extended so as 
to apply to a finite rectangular wing oscillating in a supersonic stream. 
The velocity potential and aerodynamic—' torsional-moment coefficient 
based on the linearized equations of motion for small disturbances are 
derived by means of appropriate distributions of moving sources and 
doublets. The aerodynamic— torsional-moment coefficient thus derived is 
combined with a mechanical-damping coefficient to study the effect of 
aspect ratio on the undamped torsional oscillations of a finite rec- 
tangular wing. Decreasing the aspect ratio of the wing is found to 
have a highly stabilizing effect on the undamped torsional oscillations. 
Results of some selected calculations are presented in several figures. 

It is pointed out that second— order thickness effects may be of 
significance. 


INTRODUCTION 

In theoretical studies of an oscillating wing in two-dimensional 
supersonic flow Possio noted, reference 1, that under certain conditions 
a single— degree torsional instability is possible. This instability, 
also sometimes known as a type of "single— degree flutter," was briefly 
discussed by Temple and Jahn in reference 2 and has since been further 
investigated by Garrick and Rubinow in reference 3 j by W. P. Jones in 
reference 4 , and by Cheilek and Frissel in reference 5 . 

It is pointed out in reference 3 that this single— degree flutter 
is due to the wing being negatively damped in torsion and that the 
negative damping is associated with a change in sign of the torsional- 
damping coefficient; furthermore in the two-dimensional case the 
instability may take place in a low supersonic Mach number 
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range (1 <^Mach number /2.5), at low values of the frequency parameter, 
and for axis-of— rotation locations forward of a point two-thirds of the 
chord distance from the leading edge. This phenomenon may be of 
particular importance in connection with highr-speed airplanes that are 
flown in (or through) the range of low supersonic Mach numbers at which 
the phenomenon may take place. Although in general discussions this 
phenomenon is usually associated with the word "wing" the theory applies 
as well to an aileron when the aileron is considered as a separate 
degree of freedom. Hereafter in this discussion the word "wing" may 
generally be given the broader interpretation of "wing or aileron." 

The purpose of the present paper is to extend the theoretical 
investigation of the single-degree torsional instability of a wing 
oscillating in two-dimensional supersonic flow, as presented in 
reference 3 , to a finite rectangular wing and to determine the effect 
of aspect ratio on the undamped torsional oscillations of such a wing. 

It is assumed that the negative damping phenomenon in three-dimensional 
flow, like that in two-dimensional flow> is determined mainly by low- 
order terms of a low frequency and that only the effect of the first 
power of the frequency must be considered' to get a good approximation 
of its total effect. 

In order to obtain the three-dimensional velocity potential and 
aerodynamic— damping coefficient the method suggested by Garrick and 
Rubinow in reference 6 , which is briefly discussed in subsequent 
paragraphs, is applied to a thin, flat, rectangular wing performing 
slow, sinusoidal, torsional oscillations in a supersonic stream. The 
particular wing treated is such that the Mach cones emanating from the 
foremost point of each tip do not intersect the opposite tip ahead of 
the trailing edge of the wing. 

The procedure developed herein may be readily extended to apply to 
any plan form with supersonic trailing edge as long as other edges that 
might be in the regions of mixed supersonic flow are continuously straight. 
What might be a more desirable extension would be to obtain the non- 
linear effect that thickness might have on the undamped torsional oscil- 
lations of a given plan form. 

In reference 6 , Garrick and Rubinow make use of the theory of small 
perturbations to investigate the air forces on a thin finite wing 
oscillating in a supersonic stream. For convenience the boundary— value 
problems for the velocity potential for a three-dimensional surface 
(finite wing) moving at supersonic speed are classified into two types 
and referred to as "purely supersonic" and "mixed supersonic." The 
purely supersonic boundary— value problem refers to regions of flow 
where no interaction between the flow on the upper and lower surfaces is 
present. In this case the surfaces can be treated separately and the 
boundary— value problem for each surface can be satisfied by source and 
sink distributions. The source and sink distributions for each surface 
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are known functions of the plan form and profile of the wing. Conse- 
quently, the velocity potentials in the purely supersonic region can 
always he expressed in the form of surface integrals with known 
integrands . 

The mixed— supersonic boundary— value problem refers to regions where 
interaction between the flow on the upper and lower surfaces is present. 
This problem cannot be solved by a distribution of sources alone over 
the wing surface but in general can be satisfied by use of doublet 
distributions. The manner in which the doublets are to be distributed 
(or the distribution function) depends on the camber of the wing and on 
the plan form of the region of the wing where interaction between the 
flow on the upper and lower surfaces is present. In order to find a 
required distribution function it is usually possible to make use of the 
given boundary conditions and express . the distribution function as the 
unknown function in an integral equation. 

It may be appropriate to mention that Eward, reference 7* has 
recently developed, by consideration of a source distribution over the 
entire upwash field, a time-dependent velocity potential that may be 
applied to certain edge problems. In reference 8 Harmon made use of 
this development to derive some stability derivatives for thin 
rectangular wings at supersonic speeds. The velocity potential derived 
herein for the slowly oscillating case can be shown to bear a relation- 
ship to the sum of three potentials employed by Harmon, namely those 
due to (l) constant vertical motion, (2) accelerated vertical motion, 
and (3) pitching motion. 


SIMBOLS 


disturbance— velocity potential 


V 

K 

K 

K 

x, y, z 


potential function of a moving source defined in 
equation (7) 

potential function of a moving doublet defined in 
equation (l6) 

potential function due to distribution of sources in 
region N of figure 1 

potential function due to distribution of doublets in 
region T of figure 1 

rectangular coordinates attached to wing moving in nega- 
tive x-direction 
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m 


function defining mean ordinates of any chordwise section 
of wing such as at y = y 1 as shown in figure 1 


w(x,y ,t) vertical velocity at surface of wing along chordwise 

z=0/ 


section at y = y^. ({^ 




X 


0 


t 

a 

d 

V 

c 

M 


rectangular coordinates used to represent space location 
of sources or doublets in the xy— plane 

abscissa of axis of rotation of wing (elastic axis) as 
shown in figure 1 

time 

angle of attack 

time derivative of a 

velocity of main stream 

velocity of sound 

free— stream Mach number (V/c) 


p = \/m 2 - 1 

Tgj ^ ^ >unc '* : 'i ons defined with equation ( 7 ) 




W(t - T X ) 
w(t - Tg) 


U> 

b 

h 


function used to represent space variation of source and 
doublet strengths 

functions used to represent time variation of source and 
doublet strengths 


angular frequency 
one— half cord 
one— half span 


aspect ratio 


a 
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F(|,ti), f" ( ^ j tj) functions used to denote doublet distributions 
p density 

p pressure difference, measured positive downwards, defined 

in equation (24) 

p„,'p values of p in regions N and T, respectively, of 

T figure 1 


M_ 


a 


M. 


M», 


*V 


8 , 


a 


M- 

m 

r 

( 

I 


k = ^ 
V 

V 


aerodynamic moment, defined in equation ( 25 ) 


in— phase component of aerodynamic-moment coefficient 
defined in equation ( 30 ) 

out-of— phase component or (torsional-damping) aerodynamic- 
moment coefficient, defined in equation ( 31 ) 

total— torsional-damping coefficient 


mechanical-damping coefficient; ng corresponds to the 

” 0 L 

usual logarithmic decrement 

wing density parameter ( 

\4h 2 / 


mass of wing 

radius of gyration divided by 

moment of inertia per unit length of wing about elastic 
axis 



flutter coefficient 


natural angular frequency of torsional vibration about 
elastic axis 
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ANALYSIS 


Consider a thin rectangular wing moving at a constant supersonic 
speed in a chordwise direction normal to its leading edge. In accord- 
ance with linear theory the boundary— value problem for the velocity 
potential is treated at the plane of the wing. For the portion of the 
wing between the Mach cones emanating from the foremost point of each 
tip, region N in figure 1(a), the boundary-value problem- involves only 
the purely supersonic flow; and for the portions of the wing within these 
Mach cones, regions T in figure 1(a), the boundary— value problem involves 
the mixed supersonic flow. Thus the solution to the present problem pro- 
poses the use of the two types of boundary— value problems discussed in 
the introduction. 


Boundary— value problems for velocity potential- for a rectang ular 
wing. — The differential equation for the propagation of small disturbance s 
that must be satisfied by the velocity potential in both regions T 
and N (referred to a uniformly moving coordinate system as shown in 
fig. 1 ) is (equation 4 , reference 6 ) 


0 «V* if if if 


( 1 ) 


The boundary conditions that must be satisfied by the velocity potentials 
may be stated as follows; (a) In regions T and N the flow must at 
all times be tangent to the wing surface and (b) in regions T the pressure 
must fall to zero along the wing tips and remain zero in the portion of 
the Mach cones emanating from the foremost points of the wing tips not 
occupied by the wing. (As customary with linear boundary conditions the 
effects of thickness and camber are separated. In particular, in the 
nonstationary case the linear thickness 9ffect is of no significance and 
is not considered here. The important camber effect is characterized by 
the conditions that the perturbation pressure, the perturbation velocity 
in the free— stream direction, and the perturbation— velocity potential are 
all anti symmetrical with respect to the reference plane, whereas the 
perturbation velocity normal to the reference plane is symmetrical. 
Accordingly, the boundary condition of zero pressure in the side wake may 
be stated as zero perturbation velocity in the free— stream direction. In 
view of these antisymmetric characteristics it is only necessary to derive 
the perturbation— velocity potentials for either the upper or lower surface 
of the wing. The upper surface is chosen for the derivations herein.) 

With the boundary conditions stated as in the preceding paragraph, 
the difference in the two types of boundary— value problems involved 
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is in the additional conditions stipulated in condition (h) that must 
be satisfied in region To 

The condition of tangential flow may be expressed analytically as 





, r Sz m ^m 
V 8x + dt 


( 2 ) 


where is the function defining the mean ordinates at any chordwise 

section of the wing. For the particular case of the wing performing 
small sinusoidal oscillations of maximum amplitude a Q about the span- 

wise axis x = Xq, the equation for may be written as (see fig. 1(b)) 


Z 

m 


a(x - x 0 ) - oe^t ( 1 - x Q ) 


(3) 


Substituting this expression for Z into equation (2) gives for this 
case 


w(x,t) = Ya + a(x — (4) 

Equation (4) indicates that the vertical motion of a wing moving for- 
ward and at the same time performing sinusoidal oscillations about a 
spanwise axis is equivalent to the superposition of the vertical motion 
of two similar wings moving forward, one at an instantaneous angle of 
attach and the other rotating about the spanwise axis x = Xq at an 

instantaneous rate of rotation. Since the differential equation (l) 
is linear, the velocity potential satisfying equations (l) and (4) may 
be considered as the sum of two potentials: the first corresponding 

to the first term of the right-hand side of equation (4) which 
hereinafter is denoted by and the second corresponding to the 

second term of the right-hand side of equation (4) which hereinafter 
is denoted by The symbolic forms of the velocity potentials for 

the regions N and T may therefore be written as 


^Na + ^Nd 


(5) 
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and 

= ^Ta + ^Ta 


( 6 ) 


Derivation of 0^.— Reference 6 shows that the boundary— value 

problems for unsteady motion in the purely supersonic region (region N) 
ma y be satisfied by distributions of moving sources. The potential at 
any point (x,y.,z) due to a moving source of strength w(t), varying only 
with time, located at point (t^T^O) is given in equation (7a) of refer- 
ence 6. In the present notation the expression for this potential may 
be written as 


w(t ~ T l) + w(t -T g ) 


\l{" - \)C>2 - ") 

M(x — 1 ) 

\/(n - Th ) (ip - 1) 

cp 2 

pc 

M(x - £ ) 

. \Z(i - 11K 12 - 1) 

cp 2 

Pc 


f(x - 1 ) 2 - P 2 z 2 


q 2 = y + i \/(x - I ) 2 - p 2 z 2 


( 7 ) 
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The integral form of the velocity potential at any point (x,y } z) due to 
a distribution of such sources over the | plane is 


0(x,y,z,t) 


1 

*2|3rt 


px— pz 


JO 


-1112 

w(t,n)0i dT i 

Jtii 


( 8 ) 


where W(£,i}) represents the space variation of source strength. 

The limits and are the intersections of the g rj— plane 

and the upstream-opening Mach cone with vertex at point (x,y,z). Thus 
the velocity potential at point (x,y,z) is affected only by the sources 
in the It]— plane that are within this Mach cone. (See fig. 2(a) for 
the limiting case z = 0.) 

In order to obtain for an oscillating flat wing., the time 

variation of source strength w(t) is defined as follows: 


w(t) = e lcDt 


This form for w(t) gives for the numerator of equation ( 7 ) 
\ * 

/. ^ x /, - N io)(t— ' T i) io)(t— ' To) 

v(t — ^ 2 _ ) w(t — ^ 2 ) = ® + e 


. T2+T1 

* 1 iu) «- 

= 2e la *k e 2 cos a> — — — i (9) 

2 


The space variation of source strength is defined as 


W(M) = “Va 0 


( 10 ) 
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If the expressions given in equations ( 9 ) and (10) are substituted into 
equation ( 8 ) and z is made zero, the following equation for the velocity 

potential on the upper surface of the wing, with a = a, e 1 "*, is obtained: 


iuM, 




va r p :^~ e) cos £ £n - - ,) 

P ”J 0 iu 1/(1 - la) (il2 - 1) 


dT) d! 


( 11 ) 


The integrations indicated in equation (ll) can be carried out in 
bbe form Ox a series of Bessel functions as in reference 3i however, this 
is not required here since, as is pointed out in the introduction, only the 
first-order effects of the frequency a) are necessary in the present 
discussion. The terms involving higher powers of the frequency can be 
deleted by expanding the integrand of equation (ll) into a power series ' 
in a) and dropping the terms involving this parameter to powers higher 
than the first. When this expansion is made, the indicated integration 
in equation (ll) can be carried out in closed form. There is obtained 


0Na * 



*2 l-^f(x-l) 

— - cB 2 

1 ~ \)(\ - ") 


dq d£ 



or, if the out— of— phase component 


io» is denoted by a. 



VaM 2 
2^c 


( 12 ) 


Note that even though the expression for 0 Wa in equation (12) does 

not contain the integrated effects of higher-order terms in- o>, it remains 
sinusoidal with respect to time* 



MCA TN 1895 


11 


In order to obtain 0^ it is only necessary to change the definition 
of in equation (10). In this case 


w ( £jT)) = ioxx Q ( £ - x Q ) 


( 13 ) 


Substituting the expressions from equations (9) and (13) into equation (8) 
and letting z approach zero yields 



which, after expansion to first power of oo, yields 



Equations (12) and (15) may be substituted into equation (5) to obtain the 
complete expression for 0^ at the upper surface of the wing. 

Derivation of 0 T .~ As pointed out in the introduction, distributions 

of doublets may be used to satisfy the boundary— value problem in the region of 
mixed supersonic flow (region T). The type of doublet required is that with 
axis normal to the reference plane. It will be recalled that the potential 
of such a doublet may be obtained from the, potential of a source, located 
in the reference plane, by a partial differentiation of the source 
potential with respect to the direction normal to the reference plane. 
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If such a differentiation is applied to the potential of a moving source 
in the plane (equation ( 7 ) ) , the following expression that may he shown 
to satisfy equation (l) is obtained: 



W (£ 


w(t — T l) + w(t — t 2) 

\ j ( 1 - ' ii )(’>2 - ') 


where now W( £ ,t] ) refers to doublet strength. 


(16) 


In order to obtain ^p a for the oscillating flat wing it is convenient 
to assign to the functions w(t — T i) + w(t — T 2 ) and W(|.,n) the 
expressions given in equation ( 9 ) and equation ( 10 ), respectively. 
Substituting these expressions into equation (16) gives 


02 “ 


d 

5? 


e or 


cos 


— 2Va 


03 

Pc 




VC - ii)(i 2 - ■>) 


laM 

- 2 Va 0 2 ze 


<x-0 


03 


3ln g ]/(" - ~ ^ 

p2 ° ( ,l - " ") 


0 - - ”)] 3/2 


( 17 ) 
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Expanding the terms involving a) in equation (I7) into a power series 
to the first order of cu gives 


02 « -2Va[i 2 z 2£ 


1-^1 (1-!) 


(” ' \)(\ ~ ")| 


3/2 


= -2V32^ _ fej 


Ma 


l" ” ^ ~ *)J ^ CP |( * " ^l) (^2 - ^ 


3/2 


( 18 ) 


Each of the two terms within the braces in equation (l8) can be shown to 
satisfy equation (l) to the first order of a>. Furthermore it can be seen 
that, to the first order of a>, the potential of a moving doublet may be 
considered as a superposition of two doublets of the stationary type; one 

of strength proportional to a — and the other of strength proportional 

cp 2 

to 

2 

. cP d 

The potential ^p a composed of distributions of doublets of the typos' 
in equation (l8) may be written 


0 . _ZB* fa - HasV/ / lU,n)d ' 1 

* l cp 2 y J Jq r w 03/2 

7 J - *'l)( n 2 - 1)J 


_ VzMa 
c0« 


^F 2 (|,q)dq d| 

(1 - (n 2 - 1)] 


3/2 


(19) 
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where the integration is to be extended over the wing surface S included 
in the fore cone with vertex at (x,y,z) (see fig. 2(b) for the limiting 
case z = 0) and where Fj_ and Fq are functions describing the manner in 
which the doublets are to be distributed in the region S and must be deter- 
mined in such a manner that the boundary conditions (a) and (b) will be 
satisfied. 

The normal procedure of determining Fj. and F 2 would be to formulate 

an integral equation by imposing the boundary conditions on equation (19); 
however , this procedure is tedious and difficult. It can be circumvented 
here by making use of known solutions to similar problems. The integrand 
of the first integral in equation (19) is recognized to be the same as that 
in the integral form of the velocity potential, for a similar wing at constant 
angle of attack, and the integrand of the second integral is recognized to 
be the same as that in the integral form of the velocity potential for a 
similar wing undergoing a constant rotation about its leading edge. 
Furthermore, the boundary conditions for these cases can be stated precisely 
as conditions (a) and (b), stated herein, for the oscillating wing. The 
distribution function and the velocity potential at the wing surface for 
the rectangular wing at constant angle of attack can be readily obtained 
from the known distribution of lift for this case. The distribution 
function and velocity potential at the wing surface for the rectangular 
wing undergoing a constant rotation about its leading edge can be obtained 
from the wing at constant angle of attack by linear superposition. The 
distribution functions and velocity potentials at the wing surface for 

both these cases are derived in the appendix. The value of 0 

r Ta 

(equation 19) at the wing surface can be easily deduced from these results; 
that is. 


0Ta( x >y> t ) 


— fa 

P« \ 


Max' 

cp 2 > 


l/py(x - 0y) + x sin*” 1 v [&L 

x 


VMa 

cP 3 Jt 


— — - \/ P y(x - Py) + x 2 sin 1 »/^ 

3 V x 
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2Va 

pit 



x — Py) + x sin 1 



^L_ ( x + 2py) \/py(x - Py) + 3x 2 sin 1 \/^ 

3*P 3 L V x _ 


( 20 ) 


In order to obtain 0 T< ^ it is convenient to define W(|,tj) as in 

equation (13). Substituting equations ( 9 ) and. (13) into equation (16) 
yields 



-2a(| - x 0 ) 


layi / t \ 

1 CQ3 £ ^(t) - - rj) 

ft-'#'*-') 


( 21 ) 


Carrying out the indicated differentiation and expanding to the first 
power of o> results in 


0 2 = -2p 2 za — 


■[(’• " "ON ' n j 


3/2 




= 2p 2 za< 


(n - »,) (n 2 - 1) 


3/2 


V 



( 22 ) 


The two expressions in the right side of equation (22) have the same form 
as the expressions in the right side of equation (l8). The following 
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integral expression for 0rp^ is therefore of the same form as the 
integral expression for <f ) ^ (equation 19) and the functions F^ 
and F^ can be determined by the procedure given for determining Fj_ 
and F2 : 


^Td " 


pzdxQ 


d.T) 

, Pza l 

[ / 6F 4 (|,n)d6 dq 

[(" - "i)(" 2 - ') 

3/2 Tt J 

Js j 

(1 

3/2 


The expression for at the wing surface can therefore be deduced from ^ 

the results in the appendix to be 


0Tct = 


2oxq 



x — Py) + x sin 


-1 



2d 5x - 2Py 
Pit s 


\/py(x - Py) 




■ & [ — ~ 610 * i(» - 2 ^) sin -1 yg 

(23) 

Equations (20) and (23) may be substituted into equation ( 6 ) to obtain (fa 
at the wing surface. 


Aerodynamic moment .— The local pressure difference between the upper 
and lower surfaces may be written as 



( 24 ) 
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The aerodynamic moment on the entire wing about the axis x = Xq, measured 
in the positive clockwise direction, is then (see fig. 1) 


Ma = 2 



+ 2 



(25) 


If all linear dimensions involved in equations (24-) and (25) are referred 
to the chord 2b, these equations become, in a sense, nondimensional and 
can be written 



Putting in the appropriate limits of integration (see fig. 2) results in 
the following equation for the total aerodynamic moment : 


M = 16b 3 
a 


m 


Jo 


J - * 0 )% dy ta + 16b 3 


(* 



If the nondimensional forms of and p^ are computed from 

equations (12), (15), (20), and (23) and substituted into equation (28) , 
there is obtained a complex expression for that may be written as 


Mq, = -8pb ^V^k^aA + iM^ 


(29) 


\ 
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where 



is the. in— phase component of the aerodynamic-moment coefficient, an d 

2P 3 (l _ 3X() + 3 xq 2) - 0(2 - 3*o) - - 2 - X Q^ 3Z n ~ 2 ll ( 4x Q ~ 3 ) 

2A 

(31) 

is the out— of— phase or torsional— damping aerodynamic-moment coefficient. 

The expression for in equation (31) represents the central result 

of the present investigation because the conditions of torsional stability 
or instability depend directly on the sign and magnitude of M^. It will 

be noted, in equation ( 31 ), that may undergo changes in sign as any 

one of the parameters |3, Xq, and A is varied. 

Although equation ( 31 ) has been derived for the condition that the 
intersection of the Mach cones from the foremost points of the side edges 
is off the wing, it can be shown that it remains valid as long as the Mach 
cone from one side edge does not intersect the opposite side edge. When a 
side edge is intersected by the Mach cone from the opposite edge different 
velocity potentials from those explicitly derived herein are required for 
portions of the wing behind the points where the Mach cone intersects the 
edge. The derivation of these potentials is not considered here. 

Total torsional— damping coefficient and some selected calculations of 
stability, conditions .— If mechanical damping is assumed to exist about the 
axis of rotation, for example as the hinge friction in an aileron instal- 
lation, and if this damping is converted to coefficient form as in 
references 3 and 9 and combined with the aerodynamic damping, the total 
torsional-damping coefficient may be written as 

V = m 4 (p,xo,m) + ^ a 2 (^) 2 s a 
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or 


V = M u (p >Xo ,k,A) + i nr a 2^^ 2 g^ (32) 

Negative values of M^' indicate dynamic instability that may 
correspond to single-degree torsional flutter. Positive values of M^' 

indicate stable conditions and between the stable and unstable conditions, 
that is when M^’ vanishes, a borderline state of unstable equilibrium 

separating damped and undamped torsional oscillations occurs. Since the 
term involving in equation (32) is always positive, negative values 

of M^' and the vanishing of M^' are associated with negative values 
only of M^ (equation (31)). 

With the exception of aspect ratio A, the effects of the individual 
parameters in equation (32) are presented in references 3 and k. The 
effect of aspect ratio for some particular values of the other parameters 
is briefly discussed in the following paragraphs. 

The range of values of Xq and 3^ (and M 2 ) for which M^' vanishes 

when g a is assumed to be zero and A has different selected values is 

shown plotted in figure 3» The regions inside the different curves give 
the range of values of Xq, 3 2 , and M 2 in which instability might occur. 

The curve for infinite aspect ratio agrees with that given in figure 21 of 
reference 3 and that given in figure 2 of reference h. It will be noted 
that decreasing the aspect ratio has a highly stabilizing effect. 

The dashed portions of the curves in figure 3 correspond to aspect 
ratio and Mach number combinations that cause the Mach cone from the 
foremost point of one tip to intersect the opposite tip ahead of the 
trailing edge. 

In figure 4 the quantity M^ multiplied, for convenience, by 33^k is 

plotted as a function of M 2 for the values of A selected in figure 3. 

The ordinates of these curves are proportional to the amount of negative 
aerodynamic damping available for given Mach numbers and for given values 
of the frequency parameter k. The stabilizing effect of decreasing aspect 
ratio is apparent in this figure. 

Corresponding values of the flutter-speed coefficient V/bo^ and 
aspect ratio A for which M^' vanishes for some selected values of the 

parameters Xq, go,, M 2 , p, i-, and r a 2 are plotted in figure 5. 
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The regions above these curves give the range of values in which insta- 
bilities might occur. 

It should be recalled that thickness effects have not been considered 
in the present discussion. Second— order thickness effects may be of 
considerable significance in regard to undamped torsional oscillations of 
wings of finite or infinite aspect ratio. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

• Langley Air Force Base, Va., March 28, 1949 
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APPENDIX 

VELOCITY POTENTIALS FOR TIP REGION OP RECTANGULAR 
WING IN STEADY SUPERSONIC FLOW 


It can be shown, as in reference 10, that at points in the xy— plane 
the perturbation-velocity potential in a supersonic stream due to doublets 
in the xy— plane distributed according to a distribution function f(|,T|) 
is proportional to the distribution function; that is. 


0(x,y,O) = lim 
z-jO 


Va0z 


M)dt) <3- 1 


fx - |) 2 - 0 2 (y - tj ) 2 - p 2 z 2 j 3 / 2 


±rtf(x,y) 


(Al) 


where the positive sign applies to the upper surface of the region S and 
the negative sign applies, to the lower surface; that is, the direction from 
which z 0. When the appropriate distribution function can be found, such 
doublet distributions can be used to satisfy boundary— value problems, in the 
xy— plane, for vanishingly thin airfoils at vanishingly small angles of 
attack. In the case of a rectangular wing placed in a supersonic stream 
flowing from the negative x— direction, with its leading edge normal to the 
free— stream direction, the added velocity in the free— stream direction in 
the tip region can be obtained from the expression derived for the added 
pressure difference, in this region, in reference 11. In the present 
notation. this expression for pressure difference is 


2 ,vMJ S fif sin-1 M 
8x 0 n V x 


or 


& 

8x 




(A2) 
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If 0 is eliminated from equations (Al) and (A2 ) , the following differ— 
ential equation for f (x,y) is obtained: 


* £ = 22 sin-1 J£ 
8x P 1 * V x 


(A3) 


The expression for rtf(x,y) or 0(x,y) can now be obtained by a partial 
integration, with respect to x, of equation (A3). This integration gives 


nf(x,y) = 0(x,y) = ^ 

P* 


\/py(x - 0y) + x sin 1 


(Ah) 


For the rectangular wing rotating about its leading edge the 
distribution function for the tip region can be obtained by a linear 
superposition of the function given in (Ah); that is. 


“* jf fU - £)<*£ = 



\/pyJ(x - |) - pyj + (x - |) sin ~ X 


d| 


Va 

P* 


+ X 2 sin - 1 


(A5) 


If in equation (Ah) a is replaced by a - S and if in equation (A5) 

Ma ° P 

a is replaced by — *■, the sum of the resulting equations will yield 

cp^ 

equation (20) of. the text. Similarly, if in equation (Ah) Va is 
replaced by -Xgd and if in equation (A5) Va is replaced by d, the 

sum of the resulting equations will yield equation (23) of the text. 
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(a) Purely supersonic region. 



(b) Mixed supersonic region. 

y 

Figure 2.- Sketch illustrating areas of integration for "purely supersonic" 
and "mixed supersonic" regions of flow* 
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Figure 3 .- Curves showing ranges of M 2 , and x 0 for which the aerodynamic torsional 

damping vanishes for some selected values of A • g rt se 0„ 
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Figure 4.- Continued 





Figure 4.- Continued 
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(e) x Q t 0.6. 
Figure Concluded ; 
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stable region 









Figure 5.- Continued. 














